Field- Tuned Superconductor-Insulator Transition in BaPbi_ x Bia;03 
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BaPbi-^Bi^Os is found to exhibit a field-tuned superconductor to insulator transition for Bi 
compositions 0.24 < x < 0.29. The magnetoresistance of optimally doped samples manifests a 
temperature-independent crossing point and scaling of the form p(T,H) = p c F(\H — _r7 c |T -1 ' z "), 
where H c is the field determined by the temperature-independent crossing point, and zv — 0.69 ± 
0.03. High resolution transmission electron microscopy measurements reveal a complex intergrown 
nanostructure comprising tetragonal and orthorhombic polymorphs. Data are analyzed in terms 
of both a classical effective medium theory and a field-tuned quantum phase transition, neither of 
which provides a completely satisfactory explanation for this remarkable phenomenology. 
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I. INTRODUCTION 

BaBiOa is a charge density wave (CDW) insulator 
comprising two distinct Bi sitesii - — Chemical substitu- 
tion of either K for Ba or Pb for Bi results in suppression 
of the CDW state, and for a limited range of composi- 
tions yields a superconductor, with maximum T c values 
of 12 and 30 K for BaPbi_ x Bi x 3 and Bai_ x K x Bi03 
respectively^— Given their moderately high T c values, 
these materials have attracted considerable attention in 
order to determine the superconducting pairing interac- 
tion, and factors determining the phase diagram. Signif- 
icantly, chemical substitution not only modifies the elec- 
tronic structure^ but also introduces disorder, and with 
it a whole new set of issues that can play an important 
role in both the superconducting properties and also the 
compositionally-tuned metal-insulator transition. 

Several previous studies of BaPbi_^Bi x 03 including 
x-ray diffraction^ - — transport measurements^ opti- 
cal reflectivity; 15 i 16 and thermodynamic measurements^ 
have provided important insights to the evolution of the 
crystal structure and physical properties across the phase 
diagram. The material superconducts for compositions 
0.05 < x < 0.35, but the volume fraction is always less 
than 100 %. Recent high resolution x-ray and neutron 
scattering measurements have revealed the presence of 
two polymorphs in this compositional range, one with a 
tetragonal symmetry and one orthorhombic The su- 
perconducting volume fraction scales with the relative 
proportion of the tetragonal phase, leading to the con- 
clusion that only the tetragonal phase harbors supercon- 
ductivity. However, it is not clear what factors determine 
T c in this material, and previous studies have neglected 
the profound effects of localization and enhanced inter- 
actions arising from disorder. 

In this paper we present magnetoresistance mea- 
surements of single crystal samples of superconducting 
BaPbi-zBirOs that reveal a remarkable temperature- 
independent crossing point for compositions close to 
optimal doping. To account for this phenomenology 



within an effective medium theory (used to account for 
the dimorphic nature of the material in this composi- 
tional range) requires a possibly unphysical degree of 
fine-tuning. On the other hand, the existence of a 
temperature-independent crossing point also motivates 
analysis of the data within the context of a field-tuned 
quantum phase transition (QPT). The data scales re- 
markably well for an effective dimensionality d — 2, yield- 
ing a product of critical exponents of zv = 0.69 ± 0.03. 
It is not clear why the resistivity should scale in this 
manner for this fundamentally 3-dimensional material, 
though the result is probably related to the structural 
dimorphism revealed by our high resolution TEM mea- 
surements. This surprising phenomenology exposes the 
need for a theory incorporating aspects of both percola- 
tion and quantum critical fluctuations. Our results also 
imply a key role for disorder in determining T c in this ma- 
terial, perhaps also accounting for the lower maximum T c 
of BaPbi_ K Bi x 03 relative to that of Bai-^K^BiOa. 



II. EXPERIMENTAL METHODS 

Single crystals of BaPbi-^Bi^Oa with < x < 0.29 
were grown using a self-flux technique, similar to that 
described in Ref. 0. The cation composition was de- 
termined by electron microprobe analysis. These mea- 
surements revealed a uniform composition across each 
sample, within the experimental uncertainty (± 0.02). 
T c and the superconducting volume fraction were char- 
acterized by measurements of the field-cooled and zero- 
field-cooled susceptibility using a commercial Quantum 
Design SQuID magnetometer. 

Resistivity measurements were made using a standard 
four-probe configuration. Crystals were initially cut and 
cleaved to obtain approximately rectilinear bars with an 
arbitrary orientation with respect to the crystal axes. 
Electrical contact was made using sputtered gold pads. 
To minimize uncertainty in absolute values of the resis- 
tivity due to geometric factors, at least four samples for 
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FIG. 1. (Color online) Variation of the critical temperature T c 
(left axis) and zero-field-cooled magnetic susceptibility (right 
axis, red triangles) of single crystals of BaPbi-^BirOa as a 
function of the Bi concentration x. T c values were obtained 
from measurements of resistivity (open black squares for a 
90% criteria, and filled black squares for a 50 % criteria) and 
magnetic susceptibility (black filled circles). 
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FIG. 2. (Color online) Temperature dependence of the resis- 
tivity of BaPbi-^BiajOa. Curves are labeled according to the 
Bi concentration, x. 



each composition were measured and the average value 
used to scale the curves that are shown. 

High resolution transmission electron microscopy 
(HRTEM) measurements were taken for samples with 
bismuth composition of 24%. Samples were crushed in 
liquid nitrogen cooled ethyl alcohol, and the liquid was 
allowed to warm to room temperature. The slurry was 
stirred and a small droplet was placed on a holey car- 
bon grid and dried in air. The sample was analyzed us- 
ing a FEI G2 F20 Tecnai STEM operated at 200 keV. 
Thin areas were analyzed with selected area diffraction, 
energy dispersive spectroscopy and high resolution imag- 
ing. Thin areas were aligned with the [010] zone axis 
based on indexing to the Ibmm structure (space group 
#74), showing clear lattice fringes in the HRTEM. 



III. RESULTS 

T c values obtained by magnetic susceptibility mea- 
surements were estimated using a 1% criteria (0.01 of 
the value at the base temperature of 2 K). The low- 
temperature value of the zero-field-cooled susceptibility 
(47rx(T=2K)) was used to provide an indication of the 
superconducting volume fraction. A maximum value of 
~ 25% was obtained for the optimally doped samples 
(x w0.25), this value decreasing rapidly for lower or 
higher Bi concentrations (see Fig. [TJ. This behavior is 
consistent with results previously found in polycrystalline 
samples 14 , though with a reduced maximum volume frac- 
tion, possibly due to the different conditions and cooling 
rates used to grow single crystal samples. 

Figure [2] shows curves of electrical resistivity vs. tem- 
perature at zero field for representative samples of differ- 
ent Bi compositions. T c values obtained from these mea- 
surements (using both 50 and 90% criteria), are shown 
in Fig. [T] and agree with values obtained from sus- 
ceptibility measurements. Low Bi concentrations show 
"metallic-like" resistivity, with dp/dT > 0, but very 
high absolute values of the resistivity, whereas concen- 
trations of 24% and above show "insulating-like" behav- 
ior, with dp/dT < 0. The separatrix between insulating 
and metallic behavior appears to occur very close to op- 
timal doping and before the composition where the gap 
associated with the CDW opens (x » 0.35)^ 

Magnetoresistance measurements were made for fields 
up to 14 T oriented perpendicular to the current. Curves 
of resistivity as a function of magnetic field for different 
temperatures were taken for multiple samples from each 
of the compositions studied. The data were symmetrized 
from positive and negative field sweeps, although the Hall 
component was very small for all insulating compositions. 
For metallic-like samples (x < 0.19), the resistive transi- 
tion is uniformly suppressed for increasing field or tem- 
perature, as anticipated. Results for the insulating-like 
samples (x > 0.24) are qualitatively different. Represen- 
tative data for four different compositions with progres- 
sively larger Bi concentrations are shown in Fig. [3] In 
the presence of a finite magnetic field, these samples show 
a clear crossover from a superconducting-like resistivity 
(decreasing resistivity with decreasing temperature) to 
an insulating-like one (increasing resistivity with decreas- 
ing temperature). Remarkably, the samples closest to 
optimal doping (x pa 0.24 - 0.25) show a temperature- 
independent crossing point. Plotting the crossing field 
Hcrossing for adjacent isotherms as a function of the av- 
erage temperature (T ave ) reveals the convergence to a 
"critical" magnetic field H c (Fig. [4j within the resolution 
set by the signal-to-noise ratio of individual field sweeps. 
The range of temperatures over which this crossing point 
is observed to be temperature independent varies slightly 
from sample to sample. For the examples shown in Figs. 
|2Ia) and|31[b) it is found for temperatures from approxi- 
mately 8 K down to 2 K, the lowest temperature to which 
these samples were measured (solid symbols in Fig. 2]). 
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FIG. 3. (Color online) Isotherms of resistivity vs. magnetic field for bismuth concentrations of (a) x—0.24 (2 K to 10 K in 1 K 
steps), (b) x=0.25 (2 K to 10 K in 1 K steps), (c) £=0.28 (1.5 K, 2 K, 2.5 K, 3 K, 3.5 K, 4 K, 5 K and 6 K), (d) 2=0.29 (1.5 
K to 4 K in 0.5 K steps). For the nearly optimally doped compositions x = 0.24 and 0.25 a temperature-independent crossing 
point is evident for temperatures below approximately 8 K. 




FIG. 4. (Color online) Temperature dependence of the cross- 
ing field H crosa i n g for adjacent isotherms, obtained from mag- 
netoresistance measurements. T ave is the average tempera- 
ture. Compositions that exhibit a temperature independent 
crossing point H c (x = 0.24 and 0.25) are shown by black 
solid squares and blue solid circles respectively. Open sym- 
bols show low temperature deviation from H c for some sam- 
ples measured. For higher Bi concentrations (x = 0.28 and 
0.29; triangles and stars respectively) the crossing field does 
not converge to a single value over any range of temperature. 



Other samples with the same nominal composition reveal 
a similar crossing point from 8 K down to just 3 or 4 K 



(open symbols in Fig. 



below which H cross i ng devi- 



ates from H c . In contrast, for samples with a somewhat 
higher bismuth concentration of x — 0.28 and 0.29 (Fig. 
[3fc) and (d)), the crossover from the superconducting to 
the insulating state is not demarcated by a temperature- 
independent crossing point over any range of tempera- 
tures, and Hcrossing in Fig. 2] does not converge to a 
single value. 

To gain insight to the origin of the unusual magnetore- 
sistance, HRTEM measurements were made for samples 
with a bismuth concentration of x=0.24. Representative 
results are shown in Fig. [5] The high resolution im- 
age (Fig. Efa)) reveals a well-ordered structure. A fast 
Fourier transform (FFT) of this image (Fig. [SJb)) reveals 
peaks from both tetragonal (hkl even) and orthorhombic 
(hkl even and odd, in the tetragonal notation) phases. 
By systematically masking these diffraction peaks and 
performing an inverse Fourier transform it is possible 
to recreate the spatial separation of these two distinct 
polymorphs (Fig. [5]). The result of applying a mask 
such that only the (200)t and (002)t peaks common to 
both the tetragonal and orthorhombic phases are consid- 
ered is shown in Fig. [6][a). An ordered array of planes 
of atoms is evident, indicating that the two structural 
phases are highly coherent through the crystal lattice. 
In contrast, the result of applying a mask to the (101)^ 
peaks attributed only to the Ibmm orthorhombic phase, 
shown in Fig. HJb), reveals a spatial variation due to the 
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FIG. 5. (Color online) (a) HRTEM image for a sample with Bi composition of £=0.24, seen along the [010] zone axis, (b) 
Fast Fourier transform (FFT) of the HRTEM revealing the weaker {101}t reflections (in blue squares) attributed to the Ibmm 
orthorhombic phase as well as the {200}t reflections (in red circles) from both tetragonal and orthorhombic phases. Inset 
shows an expanded view to the (200) peak, where streaks in the [102] T and [— 102] T directions can be seen. 




FIG. 6. Images reconstructed from the HRTEM after Fourier filtering via masks applied to the (a) (200)t and (002)t main 
reflections, and (b) (101)t weaker reflections of the FFT shown in Fig. [5jb). (c) Filtered enhanced image resulting from 
including both kind of reflections. 



densely intergrown nanostructure. The clusters are ap- 
proximately 5-10 nm in size, which is comparable to the 
coherence length of 8 nm calculated from H c2 (0). Finally, 
Fig. injc) shows the result of applying a mask that consid- 
ers both groups of reflections ({n00}T and {00ti}t with 
n = 2, 4, and {m0m}*r with m = 1, 2, ...), more clearly 
revealing the nanostructure. The material appears to 
comprise a network of interpenetrating orthorhombic and 
tetragonal phases. 

In addition to the diffraction peaks, weak streaks are 
also evident in the FFT shown in Fig. [SJb) running in the 
[102] T and [— 102] T directions. The streaking is weaker 
for peaks associated purely with the orthorhombic phase, 
so is tentatively ascribed to a property primarily of the 
tetragonal clusters. These features are indicative of small 
domains whose interfaces are orthogonal to the beam 
(in the plane of the sample). While quantifying these 
is very difficult, the qualitative interpretation is straight 
forward; the diffuse pattern, i.e., the streak, is orthogonal 
to the interface. In the plane of the image, the tetrag- 
onal regions therefore correspond to elongated domains 



with the long axis oriented along the [201] T and [— 201] T 
directions. These regions are evident in Fig. |5{c) as ex- 
tended regions of light contrast running approximately 
from top left to bottom right and bottom right to top 
left of the image. It is not clear, however, what the full 
3D geometry of these regions is, or their connectivity in 
the b-axis direction. 



IV. DISCUSSION 

The apparent temperature-independent crossing point 
observed for nearly optimally doped compositions is re- 
markable and requires explanation. The granular nanos- 
tructure revealed by HRTEM measurements suggests two 
initial avenues to understand the experimental results, 
based first on a classical effective medium theory, and 
second on a field-tuned quantum phase transition. We 
discuss each of these approaches below. 
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FIG. 7. (Color online) Scaling analysis of representative samples with bismuth composition of (a) x—0.24, (b) x—0.25. zv was 
initially calculated from a fit to (dp/dH)\Hc vs T (lower left inset to each panel), and refined based on regression of the entire 
data set (upper right inset) as described in the main text. Best fit values for the samples shown are 0.67 ± 0.01 and 0.72 ± 
0.01 for a;=0.24 and 0.25, respectively. 



A. Effective medium theory 

One possible scenario giving rise to critical phenom- 
ena and a possible temperature independent crossing 
point in magnetoresistance curves is percolation in a 
superconductor-non-superconductor composite. From a 
Symmetrical Effective Medium Theor y 19 ] 20 , the resultant 
resistivity p m of a composite where one of the compo- 
nents has p\ — (superconductor) and fills a fraction p 
of the space, and the other has resistivity pi and fills a 
fraction 1 — p is p rn = (l/p c )(p c — p)p2 for p < p c and 
p m = for p > p c , where p c is the percolation threshold 
(which depends on the geometry of the space). In the 
case of BaPbi_ x Bi x 03, p 2 would represent the effective 
resistivity of the material in the normal state, given that 
the superconducting part associated with p\ can become 
normal above a certain field or temperature. Therefore, 
if we assume that both p 2 and the superconducting frac- 
tion p depend on temperature and field (to account for 
the effect mentioned before, as well as effects such as vari- 
ations in T c and/or Josephson coupling between differ- 
ent superconducting regions) , we can write p m (T, H) — 
(1/Pc) \pc-p(T,H)]p 2 (T,H). With this, the condition 
to have a temperature-independent crossing point at the 
critical field H c is dp m (T,H)/dT\ H = Hc = 0, or 



dp(T,H) 



dT 



p 2 (T,H c )-[ Pc -p(T,H c )[ 



dp2(T, H) 



dT 



For this condition to be true over a range of temper- 
atures requires fine tuning of the different parameters 
involved. From Fig. [2j we see that for samples near 
optimal doping, where a temperature-independent cross- 
ing point is observed, we have that both: the temper- 
ature dependence of p is very weak, so we can possi- 
bly fulfill the condition dp 2 (T, H)/8T\ h =h c « 0, and, 
if the superconducting transition is driven by percola- 
tion, p c —p(T,H c ) << 1. Therefore, the second term 
of Eqn. [T] is expected to be small for optimally doped 
samples. Hence, to get a temperature-independent cross- 
ing point requires that the first term in Eqn. Q] be very 



small, and that the two terms are finely tuned to give 
zero difference over an appreciable range of temperatures. 
The first requirement is possibly unphysical since it is 
hard to conceive a mechanism that would give rise to 
dp{T 1 H)/dT\H=H c = 0, and the second requirement is 
at best rather unsatisfying. In particular, why should 
this fine tuning occur precisely at the composition that 
also happens to yield the maximum T c ?. 

In addition, as we show below, the resistivity for fields 
on both sides of H c scales in a very specific way. Hence, 
not only is fine-tuning required to achieve a temperature- 
independent crossing point within a classical model, but 
this same level of fine-tuning is also required for all fields 
and temperatures close to H c . As such, a purely classical 
model based on an effective medium theory appears to be 
inadequate to account for the observed phenomenology. 



B. Quantum phase transition 

The temperature-independent crossing point found for 
compositions close to optimal doping implies the possibil- 
ity of critical scaling. Indeed, quite remarkably, over the 



whole range of temperatures where H c 



shows con- 



vergence for these compositions, the magnetoresistance 
curves exhibit scaling of the form 



p(T,H)= Pc F 



\H-H C 



(2) 



where p c is temperature and field independent and 
varies from sample to sample and F(y) is a universal func- 
tion. Values of the critical exponents product zv were ob- 
tained for each sample by two methods. Initial estimates 
were obtained from a linear fit to logio((i/9/(i_ff)|^f c vs 
logioT (insets to Figs. E^a) andlTJb)). Further refine- 
ment to these values was obtained by an iterative max- 
imization (in H c and zv) of the correlation coefficient 
(R 2 ) calculated from a fitting of each part of the scaling 
function (for H > H r and H < H c ) to a generic func- 
tion (ninth order polynomial). Initially H c was varied 
for values of zv obtained by a best visual collapse of the 
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data. With the optimal value of H c , zv was then varied 
to obtain the best R 2 . The results of this procedure are 
summarized in the insets of Fig. [T](a) and [TJb) for two 
specific samples. For the upper branch (H > H c ), the 
correlation coefficient docs not show a local maximum, 
whereas for the lower branch (H < H c ), it clearly shows 
one. The zv value obtained by this method is the one 
that maximizes R 2 for the lower branch, and was con- 
firmed to give the best visual collapse of the data. For 
the specific data sets shown in Fig. [71 scaling of the full 
curves gave exponents zv — 0.67 ± 0.01 and 0.72 ± 0.01 
for x — 0.24 and 0.25 respectively. Measurements were 
made for several crystals with the same nominal com- 
positions. For all the samples analyzed, zv takes values 
between 0.65 and 0.72, and the average value obtained is 
0.69 ± 0.03. 

Magnetic-field-tuned superconductor-insulator tran- 
sitions (SITs) with temperature-independent crossing 
points have been widely studied in disordered 2D thin 
films of several different materials. A variety of critical 
exponent products have been reported for such magnetic- 
field-tuned transitions: zv — 1.3 for InO ^ 21 ' 22 and 
MoGe,^ zv = 1.5 for La2-*Sr x Cu04^ zv = 2.3 for 
highly disordered InO^^ and zv = 0.7 for a-B i 25 i 26 
and a-Nb^Sii-a;, 2 - these last two cases being similar to 
our observation for BaPbi^^Bi^Oa. In general, for a d- 
dimensional system undergoing a quantum phase tran- 
sition, the expression for resistivity inside the quantum- 
critical region obtained by finite-size scaling has the form 
p{T,y) cx l; d ~ 2 F p (L T /t;), where £ is the diverging order 
parameter correlation length, and L T is a real length 
associated with the correlation length £ T oc £ 2 in the 
imaginary time axis, introduced in the quantum-classical 
mapping of the transition by the relation it/h — 1/ksT. 
Consequently, L T oc (/i/fe^T) 1 / 2 , where z is the dynami- 
cal critical exponent, introduces the temperature depen- 
dence of the scaling function, and £ oc \y — y c \~ u , where 
v is the correlation length critical exponent, introduces 
the dependence with the tuning parameter y. With this, 
p(T,y) oc T- [d - 2 V z F(\y - y^T' 1 /^)^. This expression 
has the same form as Eqn. [2] for d = 2, and a cross- 
ing point of different isotherms at the critical value y c is 
expected only for 2D systems. Such behavior has been 
extensively corroborated for the systems mentioned be- 
fore. 

Of course the principle objection to such a scaling in 
BaPbi-^Bi^Os is that our measurements are of bulk, 
3-dimcnsional single crystals. The data could not be 
fitted to a 3D quantum phase transition (ie. with a 
temperature-dependent prefactor proportional to T -1 / 2 ) 
with physically reasonable values of z for any of the com- 
positions studied. The apparent 2D scaling found for 
optimally doped samples might be a consequence of a 
"hidden 2-dimensionality" due to the nanostructure, per- 
haps in the form of a percolating network of supercon- 
ducting sheets. Within such a scenario, subtle variation 
in the nanostructure of different samples might account 
for the low-temperature deviation from 2D scaling seen 



for some samples with x = 0.24 and 0.25 (open sym- 
bols in Fig. 2]), which presumably reflects the presence 
of a finite 3D coupling below some energy scale. The 
absence of a well-defined crossing point for samples with 
x > 0.28 would imply a deviation from 2D behavior, and 
could be related to the progressive reduction in the vol- 
ume fraction of the tetragonal polymorph with increasing 
Bi concentration^ However, to obtain a temperature- 
independent crossing point in such a scenario would re- 
quire all of the intergrowths to be aligned in a paral- 
lel fashion, so that the projection of the field onto each 
"sheet" is identical, which seems unlikely based on the 
pattern seen in Fig. [6jb). Nevertheless, the streaks seen 
in the FFT in Fig. [UJb) are suggestive of a preferred 
orientation of the polymorph intergrowths, but further 
analysis is required in order to determine if such sce- 
nario is feasible. Such a description is therefore equally 
unsatisfactory as the analysis based on the symmetrical 
effective medium theory, although the empirical obser- 
vation of a temperature-independent crossing point and 
the associated scaling remain highly suggestive of a man- 
ifestation of some kind of quantum phase transition that 
mimics the behavior of a 2d system over a wide range of 
temperatures. 



It remains to be seen whether the phenomenology of 
the field-tuned superconductor-insulator transition in op- 
timally doped BaPbi^^Bi^Os is truly driven by quantum 
critical behavior. However, such a scenario could natu- 
rally account for several specific features of this material, 
and the closely related case of Bai-^K^BiOa. Specifi- 
cally, critical scaling following the form of Eq. © im- 
plies the presence of quantum fluctuations of the super- 
conducting order parameter over a wide range of tem- 
peratures. Although quantum fluctuations have a negli- 
gibly small effect for clean superconductors with a mod- 
erate carrier density, reduction of the phase stiffness due 
to reduced superfluid density and poor screening in bad 
metals can lead to significant suppression of T c relative 
to the mean-field value i 29 ' 30 Such a scenario is consis- 
tent with the large and "insulating" -like resistivity of 
BaPbi-zBLcOs. It is also consistent with the observa- 
tion that the optimal T c occurs at the composition at 
which strong localization effects become significant and 
the resistivity begins to develop an insulating tempera- 
ture dependence. Finally, this scenario could also provide 
a natural explanation for the long-standing question as to 
why BaPbi-^Bi^Os has a lower maximum T c (ss 12 K) 
than Bai-zK^BiOa (« 30 K), since, first, substitution on 
the Bi sites (which are responsible for the bands at the 
Fermi energy 4,5 ) will effectively lead to a more disordered 
state than substitution on the Ba site, and second, the 
polymorphism in BaPbi^^Bi^Oa also introduces disorder 
and localization effects, enhancing the effect of quantum 
fluctuations in the phase of the order parameter^ 
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V. CONCLUSIONS 

In summary, we have shown that optimally doped 
BaPb 1 _ a .Bi a; 03 has a remarkable magnetoresistance, ex- 
hibiting an apparently temperature-independent cross- 
ing point over a wide range of temperatures, and scaling 
reminiscent of a 2D superconductor-insulator field-tuned 
quantum phase transition. This surprising phenomenol- 
ogy cannot be accounted for satisfactorily by either a 
classical effective medium theory, or, given the 3D na- 
ture of the material, a simple field-tuned SIT in 2 di- 
mensions. Given the complex nanostructure due to the 
dimorphic nature of BaPbi-^Bi^Oa in this compositional 
range, our observations are suggestive that a full theoret- 
ical treatment must include aspects of both percolation 



and quantum-critical scaling. 
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